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Abstract. As a result of their applications in network coding, space-time cod¬ 
ing, and coding for criss-cross errors, matrix codes have garnered significant 
attention; in various contexts, these codes have also been termed rank-metric 
codes, space-time codes over finite fields, and array codes. We focus on char¬ 
acterizing matrix codes that are both efficient (have high rate) and effective at 
error correction (have high minimum rank-distance). It is well known that the 
inherent trade-off between dimension and minimum distance for a matrix code 
is reversed for its dual code; specifically, if a matrix code has high dimension 
and low minimum distance, then its dual code will have low dimension and 
high minimum distance. With an aim towards finding codes with a perfectly 
balanced trade-off, we study self-dual matrix codes. In this work, we develop a 
framework based on double cosets of the matrix-equivalence maps to provide a 
complete classification of the equivalence classes of self-dual matrix codes, and 
we employ this method to enumerate the equivalence classes of these codes for 
small parameters. 


1. Introduction 

Codes consisting of matrices over a finite field with the rank distance have been 
employed in a variety of applications, although often under different names and 
with a different focus. Recently, they have garnered attention in the context of 
error control for network coding due to their role in the construction of lifted rank- 
metric codes [i[20]. Previously, their potential use in space-time coding has been 
investigated by Grant and Varanasi 016]; in this context, these codes are known 
as space-time codes over a finite field. Finally, these codes, restricted to square 
matrices, have also been studied by Blaum et al. [5] and Roth m in the context 
of memory chip arrays and magnetic tape recording, where it is essential to protect 
against criss-cross errors; in that work, the codes were known as array codes. In 
each of these cases, the codes of interest were those consisting of matrices over a 
finite field with the relevant metric being the rank distance: the distance between 
matrices A and B is d{A, B) := rank(A—B). To distill the primary coding principles 
from each of these contexts and make precise the types of codes we investigate here. 
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we introduce the unifying terminology matrix codes, for which we will only consider 
the rank distance. 

Given the many contexts in which matrix codes arise, it is essential that we bet¬ 
ter understand the structure and distance properties of these codes, with a focus on 
codes that are both efficient, i.e. have high rate, and effective at error correction, 
i.e. have high minimum distance. The primary focus of this work is to provide a 
framework for classifying matrix codes based on these defining properties and to 
perform such a classification for the class of matrix codes that are self-dual. We 
restrict our focus to self-dual codes because the analogue of the MacWilliams Iden¬ 
tities for matrix codes demonstrates that the inherent trade-off between dimension 
and minimum distance for a code is reversed for its dual code: if a code has high 
dimension and low minimum distance, then its dual code will have low dimension 
and high minimum distance. 

In Section!^ we begin with the notion of duality for matrix codes, highlighting 
the relationship between this and the notion of duality in the block code case. 
We then turn toward classifying self-dual matrix codes in terms of their structural 
and distance properties. This leads us to review previous work on the notion of 
equivalence for matrix codes and on the collection of linear equivalence maps for 
matrix codes. 

Next, we refine the notion of equivalence maps to apply specifically to self-dual 
matrix codes. This focus on self-dual matrix codes necessitates an additional prop¬ 
erty in our definition of equivalence, namely that an equivalence map sends a self¬ 
dual code to another self-dual code. Toward this end, in Section lTTl we characterize 
the subset of matrix-equivalence maps that commute with the dual and thus main¬ 
tain the property of self-duality. 

In Secti 0 n r 4 .ll we give a matrix-code analogue of the mass formula, which enables 
one to determine when a classification of self-dual codes for a given set of parameters 
is complete. We then outline an alternative method for enumerating the equivalence 
classes of self-dual matrix codes in Section building off an approach proposed 
in [8] to enumerate inequivalent self-dual block codes using double-cosets. Finally, 
we carry out this enumeration for matrix codes of small lengths over small finite 
fields and give the results in Section [S] 

2. Background 

As described above, matrix codes arise in a number of important contexts in¬ 
dicating the need to characterize collections of matrix codes with good distance 
properties. In particular, we seek to classify the collection of self-dual codes since 
the analogue of the MacWilliams Identities mu implies they have significant po¬ 
tential for good distance distributions. To achieve this goal, we begin with the 
notion of duality for linear matrix codes, where a matrix code C C is deemed 

linear if the collection of codewords forms a vector space over F^. 

Definition 2.1 ([5]). Let C C be a linear matrix code. The matrix dual code 
of C is given by 

C'-LMat = {F g I Tr(XF^) = 0 for all AT G C}. 

We say that C is a self-dual matrix code if C = . 

To enable a comparison of the duals of matrix codes with those of block codes, 
we must introduce a natural map translating between these types of codes. 
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Definition 2.2. For X = \xij\ € the extended row vector eorresponding to 

X is the vector 

p (^ X ^ - (3^11 G Fg 

formed by concatenating the rows of the I x m matrix X. If C C is a matrix 

code, the extended block code of C is given by 

p{C) = {p{X)\X&C}Q¥^^. 

Since Tr(Xy^) = p{X) ■ p{Y), where • denotes the standard dot product, the 
matrix dual code of C is simply the inverse image under p of the standard block 
dual code of p{C), i.e. 

^ p-i (p(C')-L) . 

Thus, there is a natural correspondence between the duals of matrix codes and the 
duals of block codes. The map p also enables a notion of a generator matrix for 
a linear matrix code: we will say that C? is a generator matrix for a linear matrix 
code C if G is a generator matrix for the corresponding extended block code p{C), 
i.e. for a\\ X & C, X = p”^(x) for some x G rowspace(G). 

It is important to note that although the dual codes coincide, the distance metrics 
for block codes are different from those for matrix codes. Throughout, we consider 
only the Hamming distance for block codes and the rank distance for matrix codes. 
Given this difference in metrics, the literature regarding the distance distributions of 
the dual codes of block codes does not apply directly. In particular, one major result 
in this area is the MacWilliams Identities for block codes, which give an explicit 
relationship between the distance distribution of a code and that of its dual code. 
Although these results do not apply to matrix codes, Delsarte as well as Grant and 
Varanasi have proven analogues of these identities for the case of matrix codes Il[5]. 
An important consequence of these results is that the inherent trade-off between 
dimension and minimum distance for a matrix code is reversed for its dual code. As 
we seek matrix codes with a balanced trade-off between dimension and minimum 
distance, we turn to the collection of self-dual matrix codes. In particular, we seek 
to classify the collection of self-dual codes based on their structure and distance 
properties. We must first review a notion of equivalence for matrix codes generally, 
and then we develop a refined notion of equivalence for self-dual matrix codes. 

2.1. Equivalence of Matrix Codes. Intuitively, two codes should be consid¬ 
ered equivalent if they share all the same properties and structure. In particular, 
equivalent codes should have the same distance distribution and the same number 
of codewords, or same dimension if the codes are linear. In the case of block codes, 
the notion of code equivalence was made more tractable by restricting the dehnition 
to declare two block codes to be equivalent if and only if there exists a linear, or 
more generally a semi-linear, invertible map between them that preserves (Ham¬ 
ming) weight. For this paper, we will restrict to considering linear equivalence only, 
but all the following results can be generalized to the case of semi-linear equivalence 
as well; the interested reader may find these results in [14]. 

It is a consequence of the MacWilliams Extension Theorem [10] that the mono¬ 
mial matrices are the only Hamming-weight preserving linear maps, where a mono¬ 
mial matrix is any matrix of the form DP for some invertible diagonal matrix D 
and permutation matrix P. Thus, the collection of linear equivalence maps for block 
codes consists of only these matrices. 
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A notion of equivalence for matrix codes was first given in |15] . Two matrix 
codes are linearly equivalent if there exists an invertible linear map between them 
that preserves rank weight of all matrices in we call such a map a linear 

matrix-equivalence map. Note that this requirement of preserving rank weight for 
all matrices is stronger than simply requiring the map to preserve rank weight for the 
matrices in a given code. In the case of Hamming weight, these two notions coincide 
as a result of the MacWilliams Extension Theorem [10]; however, no analogue of 
this extension theorem holds for rank weight [T] . 

The collection of rank-preserving linear maps on consists of compositions 

of multiplication on the left by matrices in GLi(Fq), multiplication on the right 
by matrices in GLm(Fg), and, when I = m, transposition |131115j . To further the 
analogy with linear maps on block codes, we translate the action of maps on F^^™ to 
maps on F^"*; this translation will also prove useful for the enumeration of self-dual 
matrix codes. For this purpose, we return to the map p from Definition 12.21 which 
takes a matrix to the extended row vector formed by concatenating the rows of the 
matrix. In this context, we are able to describe the linear matrix-equivalence maps 
in terms of linear maps that act on the right only, and thus situate the maps in the 
same universe. Left multiplication by L £ GL/(Fg) on A £ corresponds to 

right multiplication by on a = p{A) where Im is the mxm identity matrix 

and ® denotes the Kronecker product of matrices. Similarly, right multiplication 
by M £ GLm(Fq) on A £ F|^^™ corresponds to right multiplication by (/; ® M) 
on p(A). Finally, transposition corresponds to multiplication on the right by the 
X block matrix T = [Eji]ij whose (i, block is the mxm matrix Eji, for 
1 ^ b J < ^74, where Eji is the matrix of all zeros with a single 1 in the (j, entry. 

This change of context motivates the following extension of the definition of 
matrix-equivalence maps. 

Definition 2.3. We say that / : F^™ —F^™ is a linear I x m matrix-equivalence 
map on extended row veetors if there is some linear matrix-equivalence map g : 
fixm -pixm that, for all A £ F^^™, f{p{A)) = p(g{A)), i.e. f = po g o p~^. 
We denote the collection of linear I x m matrix-equivalence maps on extended row 
vectors by Equivyed^g^"^)- 

From the previous commentary and the fact that XIi ® is the identity map 

on Fg™ for any A £ F*, it is easy to show that the collection of matrix-equivalence 
maps on extended row vectors has the group structure given in Proposition 12.41 
below, and so we omit the proof here; a full proof is available in [14] for the interested 
reader. 


Proposition 2.4. The group Equivlinear I x m matrix-equivalence 
maps on extended row vectors in F^"* satisfies 


EquiVVec 


{TdL^ ®M) I i £ {0,1}, L, M £ GLi(F,)} if I = m 
I L £ GL/(F,), M£GL^(F,)} if I d m 


where T = [Ejdj is the matrix for transposition. The group structure is given hy 


EqmvvedK''^) = 


Z 2 X {GLi (F,) X GU (F,)) /N if I = m 
(GLz(F,) xGL„(F,))/iV if I dm 


where N = {(A/;, A-iJ^) | A £ F*} < GU(¥g) x GL„(F,). 
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Remark 2.5. By Proposition 12.41 each linear matrix-equivalence map on extended 
row vectors has the form r®(L^ (g) M) although there are many maps that have 
the same action as T'^[L^ ® M) since T^{L^ ® M) produces the same effect as 
T''{L^ (g M){XIi (g for any A € F*. For simplicity, we will use the notation 

(g M)] to denote the equivalence class of maps corresponding to (g 

M).((A/z®A-U^) I A G Fp. 

3. Method for Enumerating Self-Dual Matrix Codes 

The goal of this paper is to enumerate all the linearly-inequivalent q-ary self-dual 
I X m matrix codes for small values of g, I, and m. The current definition of code 
equivalence is too weak for this purpose, however, because it is possible for self¬ 
dual codes to be considered equivalent to codes that are not self-dual. Section 13.11 
below gives a sufficient condition that may be added to the definition of equivalence 
to ensure that the equivalence classes of self-dual codes will consist solely of self¬ 
dual codes. Also in that section is a characterization of the matrix-equivalence maps 
that satisfy this additional condition together with a review of the block-equivalence 
maps that achieve this condition. 

3.1. Equivalence Maps that Commute with the Dual. Recall from Section 
o that we have required equivalence maps to be invertible, linear, and weight¬ 
preserving. Motivated by [181 p. 185], we now will also require that an equivalence 
map / satisfies /(C-*-) = (/(C'))-*- for every linear code C, since this property of 
commuting with the dual is precisely what is needed to ensure self-dual codes are 
mapped to other self-dual codes. 

The subset of linear block-equivalence maps that satisfy this property of com¬ 
muting with the dual, i.e. preserving orthogonality, will be known as linear block- 
equivalence maps for self-dual codes and denoted by LEquiv^g^i^{¥^). Recall from 
Section o that the set of general linear block-equivalence maps is Mon„(Fq), the 
group of monomial matrices, i.e. the group of invertible matrices of the form DP 
for some diagonal matrix D and permutation matrix P. We are interested in the 
following subgroup of Mon„(Fq): 

Definition 3.1. The scalar monomial group is the subgroup Mn{¥q) of Mon„(Fq) 
consisting of matrices of the form DP where P is a permutation matrix and D is a 
diagonal matrix whose diagonal entries are ia for some a G F*. 

MacWilliams shows m that the matrices in AI„(Fg) are the only monomial 
matrices that commute with the dual, and so Equiv^^^i^{¥q) = Al„(Fq). 

In the matrix code setting, we have: 

Definition 3.2. An invertible map / : —>• is a linear matrix-equivalence 

map for self-dual codes if / is F^-linear, preserves rank weight, and has the property 
that for all linear codes C C . The collection of linear 

matrix-equivalence maps for self-dual codes is denoted and we say 

that two self-dual matrix codes C, C C F^^'" are linearly matrix-equivalent if there 
exists a linear matrix-equivalence map for self-dual codes / such that C = f{C). 
Finally, the collection of linear matrix-equivalence maps on extended row vectors 
for self-dual codes is denoted Equiv^^iW^f^'^) and it consists of all maps / for which 
there exists g G Equiv^^^{¥^^™') such that, for all A G F^^™, f{p{A)) = p{g{A)), 
i.e. / = pogop-^. 
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We will use the orthogonal similitudes group to characterize the matrix-equivalence 
maps on extended row vectors that commute with the dual: 

Definition 3.3 ([22]). The orthogonal similitudes grow^ also called the general 
orthogonal group mi p. 136], is the collection of matrices 

GO„(Fg) = {A e GL„(Fg) I = A/„ for some A € F*} 

For A G GO„(Fq) with AA^ = XIn, we call A the similitude charaeter of A. The 
subgroup 0„(Fq) of GO„(Fq) consisting of matrices with similitude character 1 is 
the orthogonal group. 

Remark 3.4. The orthogonal similitudes group may also be defined more generally 
on any vector space V with a symmetric bilinear form (•, • ) (THiiiiim]. In this 
case, the orthogonal similitudes group is the collection of matrices that preserve that 
bilinear form up to a scalar, again called the similitude character. The orthogonal 
group is the subgroup consisting of matrices with similitude character 1. Note 
that the definition of the orthogonal similitudes group given in Definition 13.31 is the 
special case of this definition for D = F^ with the standard dot product. In classical 
group theory, a different symmetric bilinear form (or in the case of characteristic 
2, a quadratic form) is often used when defining the orthogonal group [T] p. 39, 
113]. When n is even there are two equivalence classes of symmetric bilinear forms, 
which give rise to non-isomorphic groups; we show in the appendix which group is 
conjugate to the group that preserves the standard dot product that we consider 
here. 


We now turn to a characterization of the matrix-equivalence maps on extended 
row vectors that commute with the dual. 


Proposition 3.5. The group Equivy^J^^^^) of linear I x m matrix-equivalence 
maps for self-dual codes on extended row vectors satisfies 


i?i7m4£(F'>^™) ^ I 

where N = {(A/i,A“^/m) | A 


Z2 K (GOi(F,) X GO0q))/N 
(GO,(F,) X GO^(F,))/iV 

e F*} < GO,(F,) X GO™(F,). 


if I = m 
ifly^m, 


Proof. By Proposition [HU] we have 

0 M) I * e {0,1}, T, M e GL^F,)}, if / = m 
^ vec ) \ {L^^M \Lg GL,(F,), M g GL™(F,)}, iil^m. 

Thus, we need only determine which maps in this subgroup commute with the 
matrix dual. Note that these maps act on p{C) for any C C go we need to 

express the matrix dual code in terms of p{C). By the commentary following 

Definition l2.21 we see that p(C-*““'^*) = (p(C))^ where T denotes the standard block 
code dual. Thus, it suffices to determine which maps in Equivyec (F^^'") commute 
with the standard block code dual. We will handle I = m and I m simultaneously 
by writing Flgmuvec (F^^”") = {T*(L^ G> M) \ i G {0,1}, L G GLi(F,j), M G 
GLm(Fg)} where i = 0 ii I ^ m. 

Let / G Equivy^^{¥'-^^'^). Then since Equivy^^{¥lj^"^) C iJq-«*uvec(Fq^™), there 
exists some L G GL;(Fg) and M G GLm(Fg) such that / = T'{L ® M), and / 


*There does not seem to be standard notation for this group, and so we have chosen to follow 
|21l I22| . The documentation for Magma refers to this group as the conformal orthogonal group 
denoted by COn(Fq). 
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satisfies f{{p{C))-^) = if{p{C)))^ for every linear matrix code C C Since 

U is a linear matrix code, there is a generator matrix Gc S such that 

p{C) = {a.Gc I a G F^} where k = dimF,(C'). Then (p(C'))-*- has Gc as its parity- 
check matrix, i.e. {p{G))^ = {b G F^™ | Gcb^ = 0}, and so 

fiipiC))^) = {br (lT ® M) I Gch^ = 0}. 

Since GcT'^{L^ ® M) is a generator matrix for f{p{G)), we have 

(/(P(C)))^ = {b I GcT\L^ ® M)h^ = 0}. 

Since /((p(G))-L) = {f{p[G)))^, each bT*(LT 0 M) G f{{p{C))^) with Gch'^ = 0 
must satisfy GcT^{L^ <8) M)(bT®(L^ (8) M))^ = 0. Hence, Gc and GcT''{L^ ® 
M){T'^{L^®M))^ must have the same nullspace for every possible generator matrix 
Gc, and so T^(L^ ® M){T^(L^ ® Af))^ = \Iim for some A G F* since any column 
operations other than global multiplication would change the nullspace of at least 
one generator matrix Gc- Thus, 

= T%L^ ®M){T\L^ ®M))^ = 

Using the facts that T is symmetric and = Iim, we have 

(Ml 0 MM^) = = Xllra- 

Since Xlim = A(J; 0 M) and decomposition into Kronecker products is unique up 
to scalars, we have that L^L = Xili and MM^ = X 2 lm where A = A 1 A 2 ^ 0. 

Thus, L G GO/(Fq) and M G GOm(IP'g), and so 

- / MHM0M) I *G{0,1}, L,MgGO,(F,)}, ifZ = m 
^ ) \ {L^ (^m\L& GOi(F,), M G GO^(F,)}, if Z ^ m. 

Finally, using the group structure of the Kronecker product and the semi-direct 
product structure for the action of the transpose matrix, we obtain the desired 
isomorphism. 

□ 

Remark 3.6. From the proof of Proposition 13.51 we see that the linear matrix- 
equivalence maps for self-dual codes correspond to the maps in (T) x {GOi{¥q) 0 GOm(]Fg)), 
a subgroup of GOim(®'g)- Meanwhile, the linear block-equivalence maps for self-dual 
codes are the scalar monomial matrices M-im{¥q), also a subgroup of GOim(IFg)- 
These two groups do not coincide in general. For example, over odd characteristic, 
when I = m = 2, the block matrix 


1 

-1 

0 

0 ‘ 




1 

1 

0 

0 

= IiG 

■ 1 

-1 ■ 

0 

0 

1 

-1 

1 

1 

0 

0 

1 

1 





is an element of GOi(Fq)0GOm(Fq) and thus also an element of GO;m(Fq); however, 
this matrix is not an element of Mim(¥q). A similar trick can be used to show the 
distinction between these groups for any other value oi I > 1. Additionally, there 
are matrices in A4imi¥q) that cannot be written as a Kronecker product, and thus, 
are not elements of GOi(¥q) 0 GOm(]Fg) or of (T) k {GOi(¥q) 0 GOm(IFg))- 
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4. Characterizing Linear Equivalence Classes 
OF Self-Dual Matrix Codes 

We now investigate the linear-equivalence classes of self-dual matrix codes of 
small parameters. A similar such analysis of self-dual block codes was an important 
achievement in classical coding theory. One practical motivation for this analysis 
arises from the MacWilliams Identities, which show that the relationship between 
dimension and minimum distance for a code is reversed for its dual code. As a result 
of this reversed relationship, many self-dual codes attain a perfectly balanced trade¬ 
off between dimension and minimum distance, making some self-dual codes the best 
codes known for given parameters m- This analysis continued for suboptimal self¬ 
dual codes as well, however, because of the nice mathematical properties of self-dual 
codes, which have produced applications in groups, lattices, and designs [12]. In the 
sequel, we follow roadmaps from the study of self-dual block codes to enumerate 
representatives of the equivalence classes of self-dual matrix codes. 

Section 14.11 focuses on one key tool that was used to enumerate self-dual block 
codes, namely the mass formula. This formula was essential for allowing researchers 
to determine if all the inequivalent self-dual block codes they had found through 
various techniques actually gave a complete enumeration of the equivalence classes. 
We give an analogous mass formula for self-dual matrix codes that is virtually iden¬ 
tical to that from the block code case, taking into account the different equivalence 
maps for matrix codes. These formulas align so closely because they are both direct 
applications of the Orbit-Stabilizer theorem and because there is a one-to-one corre¬ 
spondence between self-dual block codes and self-dual matrix codes as demonstrated 
in Section HD 

If one has already generated a number of constructions for self-dual codes and 
simply needs to confirm that all equivalence classes have been accounted for, then 
the mass formula is particularly valuable; however, this formula does not actually 
give a method for producing/finding these self-dual codes. In Section IL^ we turn 
to an alternative method for enumerating inequivalent self-dual codes that uses 
double cosets. This method will prove more useful here because it will not require 
the development of multiple new constructions of self-dual codes; we will show that 
we can construct a single canonical self-dual code for each set of parameters and 
this will suffice to produce a list of all inequivalent self-dual matrix codes. 


4.1. Mass Formula for Self-Dual Matrix Codes. As mentioned above, the 
mass formula for self-dual block codes is used to determine whether a given set of 
linearly-inequivalent self-dual codes is a complete representation of all the equiva¬ 
lence classes of such codes. This formula comes directly from the Orbit-Stabilizer 
theorem of abstract algebra, which can be found in any introductory graduate text, 
e.g. lU Proposition 2, p. II4]. For completeness of the duality theory, we give an 
analogue of the mass formula for self-dual matrix codes, although we follow a differ¬ 
ent method, described in Section 021 to actually enumerate the equivalence classes 
of self-dual matrix codes. 

We begin by observing that self-dual matrix codes are in bijective correspondence 
with self-dual block codes of appropriate dimensions. 

Lemma 4.1. Let C C he a matrix code with corresponding block code p{C) C 

Then C is a self-dual matrix code if and only if p{C) is a self-dual block code. 
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Proof. This is an immediate consequence of the commentary following Definition 
12.21 which showed that = (p(C'))-*-. 

□ 


Corollary 4.2. The number of [I x m, ^]q self-dual matrix codes is 

b n (9^+1) 

i=l 

, ^ ,71 1 r 77 ,7 ( 0 if q = 3 (mod 4) and 41 Zm 

where for even q, we set b = 1, and tor odd q, we set b = < „ . 

7 j 2 otherwise 

Proof. It was shown in [181 P- 184] that the formula above counts the number of 
self-dual block codes, and so by Lemma 14.11 the formula counts the number of 
self-dual matrix codes as well. □ 


Next we observe that the group of equivalence maps on self-dual codes acts 
on the collection of self-dual codes, and the orbits under this action are precisely 
the equivalence classes of self-dual codes. The stabilizers under the action are the 
automorphism groups of the codes restricted to the set of linear matrix-equivalence 
maps for self-dual codes. We denote such an automorphism group for a self-dual 
code C by With this terminology in place, we may now state the mass 

formula analogue for self-dual matrix codes. 

Theorem 4.3 (Mass Formula for Matrix Codes). Let {Ci C be repre¬ 

sentatives of the distinct equivalence classes of linearly matrix-equivalent self-dual 
matrix codes of size I x m over Fg. Then 


hn — i 

b n (9^+1) 


E 

linearly matrix-inequivalent Ci 


\EqmvffJ¥^^x^)\ 

\Autffat(.C.)\ 


where for even q, we set b 


1 , and for oddq, we setb = 


0 if q = 3 (mod 4) and 4 -j" Zm 
2 otherwise 


Proof. The space of self-dual matrix codes is partitioned into equivalence classes, 
or orbits, by the action of the linear matrix equivalence maps. The number of codes 
is the sum of the sizes of the orbits. By the Orbit-Stabilizer theorem, see e.g. [H 


Proposition 2, p. 114], the size of the orbit of a code C equals 
The result follows. 


\Eqmv^Z(¥‘>^"^)\ 




□ 


As noted above, the mass formula plays a crucial role in the enumeration of 
linearly inequivalent self-dual block codes. The mass formula is only useful, however, 
if one has developed a number of different constructions of self-dual codes. While we 
could take advantage of known constructions from the block code literature, we take 
a slightly different approach here that enables us to exploit a single construction 
for a canonical self-dual matrix code for each set of parameters. This alternative 
approach is described in the next section. 
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4.2. Double-Coset Characterization of Inequivalent Self-Dual Matrix 
Codes. In this section, we extend the work of Janusz [S] to recast the problem 
of enumerating linearly inequivalent matrix codes as a problem of enumerating 
double-coset representatives within a group that acts transitively on the collection 
of self-dual codes. While the problem of enumerating double-coset representatives 
is also intractable for large groups, it has proven feasible in the small cases we will 
examine here whereas a brute force search is not possible for many of these cases. 
This framework also provides a nice bridge for relating the equivalence classes of 
self-dual block codes with those of self-dual matrix codes as well as relating the 
automorphism groups of the codes in each of these contexts. Toward this end, we 
introduce a number of theorems and terminology given in Janusz’s paper for the 
case of self-dual block codes over the field F 2 , and we extend these ideas to work 
for both block and matrix codes over arbitrary finite fields. 

Janusz shows [8] that 0„(F2) acts transitively on the collection of self-dual block 
codes over F 2 . In Theorem 14.41 we prove that GO„(Fq) acts transitively on the 
collection of self-dual block and matrix codes over arbitrary fields. We include 
the proof here because of this modification, but also because the following proof is 
dramatically simpler than the proof given in [5] . 

Theorem 4.4 (Extended from [8]). The group GO„(Fq) of orthogonal matrices 
acts transitively on the set of self-dual codes of length n over F^. That is: 

1- If C C F^ is self-dual and A € GO„(Fq), then CA is self-dual, and 
2. If C, C C F^ are self-dual, then there is some A £ GO„(Fg) such that 
CA = C. 


Proof. Let G be a generator matrix for the self-dual code C. Then GG^ = Ok, 
where k = ^ is the dimension of the code. Let A £ GO„(Fg), i.e. AA^ = XIn 
for some A £ F*. Then GA is a generator matrix for CA. Since GA(GA)^ = 
GAA^G^ = GXInG^ = XGG^ = XOk — Ok, we see that CA is self-dual as well, 
and so (1) holds. 

Let G, G be generator matrices for C and C respectively that are in reduced 
row echelon form. Since up to column permutation, every code has a systematic 
generator matrix, we see that G and G have the form \Ik\M]P and [Ik\M]P, for 
some matrices M, M £ F^^^ and some permutation matrices P,PG On(Fq), where 
k = ^. Since C and C are self-dual, = MM^ = —Ik- Then GA = G for 


A = P-^ 


'Ik 

Ofc 


Since the permutation matrices satisfy PP 
A £ 0„(Fg) and C = CA, so (2) holds. 


Ok 

-M^M 

T _ J 

— -‘n 7 


p. 

we have that AA"'' = J„. Thus, 

□ 


To make use of this transitive action in our enumeration of self-dual codes, we 
must begin with some specified canonical self-dual code on which we will act. To 
develop such a canonical self-dual code when g = 3 (mod 4) , we will need a straight¬ 
forward number theoretic lemma: 


Lemma 4.5. [211 Lemma 11.1, p. 138] If q = is a power of an odd prime, then 
for any a £ F^, there exist elements a, 6 £ F^ such that -I- 6^ = a. 

In particular, Lemma 14.51 guarantees the existence of elements a, 6 £ F^ such that 
a'^ -\-b^ = —1 when g is a power of an odd prime. It is clear that for a self-dual code 
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to exist in F”, n must be even. Furthermore, Pless shows in [16] that for a self-dual 
code to exist in when q = 3 (mod 4), n must be a multiple of 4. Abiding by 
these constraints on n, we define a canonical self-dual code Ccanon in F^ for every 
value of n that admits a self-dual code. 


Definition 4.6. Let n be even and k = The canonical self-dual code Ccanon 
F” is the code generated by one of the following generator matrices as dictated by 
the value of q mod 4: 


[Ik I Ik] 
[ik I nAj 


r 

2 

0 

bik 

2 

cl k 

2 

0 

Ik 

2 

cl k 

2 

-bIk 

2 J 


where a G Fg satisfies = —1 when 
6^ -t- c^ = —1 when q = 3 (mod 4). 


if 2 I (?; 

if (7 = 1 (mod 4); 

if (7 = 3 (mod 4), 4 | n, 
q = 1 (mod 4), and b, c G F^ satisfy 


It is easy to check that for each value of q and n, the rows of the generator 
matrices given in Definition 14.61 are pairwise orthogonal and linearly independent. 
Note that different values of a, b, and c satisfying the prescribed conditions will give 
distinct, possibly inequivalent, canonical self-dual codes; however, the enumeration 
method outlined below is independent of these values, and so we may choose any 
values for these parameters to successfully carry out the enumeration. 

We now characterize the set of elements in GL„(Fq) that fix the canonical self¬ 
dual code. We use this characterization to determine which elements of GO„(Fq) 
fix that code. This is a key step to determining all the linearly inequivalent self-dual 
codes. 


Definition 4.7. Let C C F” be a self-dual code. The linear stabilizer Stab{C) of 
C is given by 

Stab{C) = {M G GL„(F,) | CM = C}. 

The orthogonal linear stabilizer GOStab[C) of C is given by 
GOStab{C) = Stab{C) n GO„(F,). 

Remark 4.8. For C = p{C) for some C G F^^™, neither Stab{C) nor GOStab{C) 
coincides with the traditional automorphism group (one extended row vectors) of 
the code C since that group is restricted to only contain equivalence maps that fix the 
code, i.e. elements of the form (T) x GOi(Fg) (8) GOm(®’ 5 ) that fix C. However, both 
Stab{C) and GOStab{C) contain the traditional automorphism group on extended 
row vectors of the code C. 


Proposition 14.91 below precisely characterizes the elements of Stab{Ccanon) for 
the canonical self-dual code Ccanon- This is an extension of Theorem 11 in |8], 
which characterizes the matrices in GL„(F 2 ) that stabilize a given binary self-dual 
code. The proof is similar to that of jS], and so we omit it here. 


Proposition 4.9 (Extended from [8|). Let Ccanon Q If'g be the canonical self-dual 
code of dimension k = ^ with generator matrix [Ik \ M] for appropriate M. Let 


Ik M 
Ik Ok 


Then the linear stabilizer of Ccanon is given by 


StabiCcanon) = { 


A 0 
B C 


S 


A, C G GLfc(Fg) and B G F^^'' J . 
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We can now precisely enumerate which elements of GOn(Fq) will map Ccanon to 
linearly block-inequivalent or linearly matrix-inequivalent self-dual codes. Theorem 
14.121 shows that the linear maps corresponding to linearly inequivalent codes are 
the representatives of distinct double cosets within the group GO„(Fq). This result 
is analogous to that found over F 2 in [8]. Before proving this theorem, we briefly 
review the group theory topic of double cosets. 

Definition 4.10. [H p. 117] Let H,K <G and let g € G. The H-K double coset 
of g is given by 

HgK = {hgk \ h G H, k G K}. 

Remark 4.11. As with left and right cosets, the collection of H-K double cosets 
partition the group G. To see this, we recast the H-K double coset of g as either 
a union over h G H oi left cosets of the form hgK or alternatively, as a union over 
k G K oi right cosets of the form Hgk. One important difference between double 
cosets and left or right cosets is that double cosets need not all have the same size, 
i.e. \HgK\ need not equal \HgK\ for g,g G G. 

Recall that GO„(Fg) acts transitively on the collection of self-dual codes, and so 
every self-dual code can be written in the form C = CcanonA for some A G GO„(Fq). 
Theorem 14.121 below demonstrates that self-dual codes CcanonA and CcanonR are 
linearly inequivalent precisely when A and B lie in distinct double cosets dictated 
by the appropriate group of equivalence maps. 

Theorem 4.12 (Extended from [8]). Let n = Im for some positive integers I and 
m, and let p : F^^™ —>• F^ be as in Definition \2.A Let GOStabiCcanon) be the 
orthogonal linear stabilizer as in Definition let A4„(Fq) be the scalar mono¬ 
mial matrices as in Definition \3.1[ and let Equiv^i.(¥y^"^) be the linear matrix- 
equivalence maps for self-dual codes as in Definition \S.A . Then for any self-dual 
codes G = GcanonA C F^ and C = CcanonB C F^ with A,Bg GO„(Fq), we have 

1. C is linearly block-equivalent to G if and only if A and B are in the same 
GOStab{CCanon)-A4n(^q) double coset o/GO„(Fq). 

2. p~^(C) is linearly matrix-equivalent to p~^(C) if and only if A and B are in 
the same GOStab{Ccanon)-Equiv^i.(¥’‘^^^) double coset o/GO„(Fq). 

Proof. By definition, the self-dual codes G and G are linearly block-equivalent if and 
only if there exists a map / G Equiv^^^/^fFg) such that C = f{C) or equivalently 
if and only if there exists sen M G Ainf^q) such that C = CM, i.e. 

GcanonA = GcanonBM. 

This equality holds if and only if Gcanon = CcanonBMA~^, which occurs precisely 
when BMA~^ hxes Gcanon, i.e. when 

BMA-^ G GOStab{Ccanon). 

There exists an A € GOStab{Ccanon) such that BMA~^ = A if and only if B = 
NAM if and only if 

B G GO>S'ta6(Gcanon)AA4n(Fq), 

i.e., A and B are in the same GOStab(Ccanon)-A4n(Fq) double coset of GO„(Fq). 

By Definition 13.21 p~^(C) is linearly matrix-equivalent to p~^(C) if and only if 
there exists an / G (F^^™) such that G = /(G). By ProDOsition l3.51 such 
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an / exists if and only if there exists an M £ Equiv^^{¥^g^^) such that C = CM, 
and so replacing M.n{¥q) with in the previous argument gives the 

analogous result for matrix-equivalence. □ 

We have reduced the problem of enumerating linearly inequivalent self-dual codes 
to the problem of determining double-coset representatives. Janusz notes that the 
problem of finding double-coset representatives in large groups is also intractable, 
so it is not clear that this reduction has significantly aided in the classification of 
inequivalent self-dual matrix codes. We find, however, that this method is more 
computationally feasible than brute force techniques, and so we are at least able to 
employ this method to enumerate the inequivalent self-dual matrix codes for small 
values of q, I, and m. The results of this enumeration are given in Section [5] 

5. Enumeration of Linearly Matrix-Inequivalent 
Self-Dual Matrix Codes 

In this section, we provide a complete enumeration of the linearly inequivalent 
self-dual matrix codes for small parameters. To derive this enumeration, we imple¬ 
mented the method outlined in Section [4.21 in the computer algebra system Magma 
using a built-in function for enumerating representatives of double cosets of matrix 
groups. These representatives give maps from the canonical self-dual code to a set 
of linearly inequivalent self-dual matrix codes. In Tables [T][6] below, we provide 
generator matrices for the linearly inequivalent self-dual codes that are the images 
of the canonical self-dual code under these maps, which by Theorem 14.121 gives a 
complete enumeration of the linear equivalence classes of self-dual matrix codes. In 
addition to the equivalence class representatives, we also provide the rank-distance 
distribution, or equivalently, since the codes are linear, the rank-weight distribution 
so that we may evaluate the optimality of these codes. We provide these descrip¬ 
tions for every equivalence class of self-dual matrix codes over F 2 ,F 3 ,F 4 , and F 5 
that was computationally feasible in Magma via the double-coset method. 

We first focus on an example of the enumeration of the linearly inequivalent self¬ 
dual matrix codes. Recall from Lemma ITT] that a self-dual matrix code C C F^^™ 
gives rise to a self-dual block code p{C) C F^'” via the map p. Furthermore, by 
definition, any generator matrix for the matrix code C is also a generator matrix for 
the block code p{C). Thus, it is natural to compare the linear-matrix-equivalence 
classes of self-dual matrix codes to the linear-block-equivalence classes of self-dual 
block codes of the appropriate parameters. Here, we consider the example of lin¬ 
early inequivalent self-dual matrix codes in Fg^^ and compare these to the linearly 
inequivalent self-dual block codes in F|. In this context, we explicitly compare the 
matrix-equivalence classes to the block-equivalence classes in order to gain further 
insight into the relationship between block- and matrix-equivalence. This example 
highlights some key differences in these notions of equivalence, which will prove 
valuable for comparing the enumerations in the block- and matrix-code case more 
generally. 

Example 5.1. Set q = 3, I = 2, m = 4. Following the method outlined above. 
Magma produces 13 linearly inequivalent self-dual matrix codes in Fg^^; the gener¬ 
ator matrices for these codes are given in Tabled In contrast, Mallows, Pless and 
Sloane show that there is only one linear equivalence class of self-dual block codes 
in Fi m- Thus, the 13 different self-dual block codes that arise from the genera¬ 
tor matrices in Table [2] are all linearly block-equivalent. To gain some intuition as 
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to why this is true, we focus on the first two linearly-inequivalent matrix codes in 
Tableland examine them both as matrix codes and as block codes. 

Let Cl C be the matrix code generated by Gi in Tabled Similarly, let 

C 2 C be the matrix code generated by G 2 in Table El Since these matrix codes 
have different rank-weight distributions, there cannot be a rank-weight-preserving 
map from one code to the other, and so these codes cannot be matrix equivalent. 

We can check that both p{Ci) and p{C 2 ) have Hamming-weight distributions 
with Aq = 1, ^3 = 16, Aq = 64 and the remaining Ai equal to zero, where Ai 
denotes the number of codewords of weight i. Thus, it is possible that there exists 
a Hamming-weight-preserving map between the codes. In fact, one can check via 
Magma, for example, that the matrix M given by 


M = 


-1 

0 

0 

0 

0 

0 

0 

.0 


0 

0 

0 

2 

0 

0 

0 

0 


0 

1 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

2 


000 

000 

002 

000 

200 

000 

020 

000 


0- 

0 

0 

0 

0 

1 

0 

0. 


eGLslFa) 


maps p{Ci) to p{C 2 )^ and since M is a monomial matrix of Is and 2s, i.e. Is and 
— Is over F 3 , M is a linear block-equivalence map. 

Intuitively, p{Ci) and p(C 2 ) are linearly block-equivalent because it is possible 
to permute and scale the entries of one code to obtain the other, and the actions 
of permuting and scaling entries do not affect the Hamming weight of a vector; 
however, these same actions will affect the rank weight of the matrix obtained from 
the extended row vector, and so the corresponding matrix codes Ci and C 2 are not 
equivalent. 

In general, when we restrict to linear equivalence maps for self-dual codes, there 
are significantly more block-equivalence maps for self-dual codes than there are 
matrix-equivalence maps for self-dual codes. In this example, applying a simple 
counting argument^ and a formula for the size of On(Fg) [ 11 ] to the groups of 
equivalence maps described in Definition 13.11 and Proposition 13.51 we see that 


lEquiv^ZckC^Dl = 10,321,920 > 18,432 = 4)| 


demonstrating that linear block-equivalence for self-dual codes is much broader than 
linear matrix-equivalence for self-dual codes. Thus, we would expect many more 
equivalence classes for self-dual matrix codes than for self-dual block codes. 

It is worth noting that while there are significantly more linear block-equivalence 
maps for self-dual codes than there are linear matrix-equivalence maps for self-dual 
codes, this phenomenon does not hold true for general linear block-equivalence and 
linear matrix-equivalence maps. If we do not restrict to those maps that commute 
with the dual, we see that we do not gain any linear block-equivalence maps because 
all monomial matrices over F 3 consist of only ± 1 , but we have drastically more linear 
matrix-equivalence maps: 

\EquivBiock{Vl)\ = 10,321,920 < 582,266,880= |£;gmz;Mat(F^^^)|. 


Thus, in the general case, we would expect to see far fewer linear matrix-equivalence 
classes than linear block-equivalence classes, even though this behavior is reversed 
in the case of self-dual codes. 


t Formulas derived from such a counting argument are available in M for the interested reader. 
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Another important question to consider is whether any self-dual codes are optimal 
with respect to their minimum distance. The analogue of the Singleton bound for 
matrix codes [3] guarantees that for any I x m matrix code, we must have 

k < min{Z, m}(max{Z, m} — d -|- 1) 

where k is the dimension of the code and d is its minimum distance. Using that 
I < m and k = this simplifies to d < ^ -I- 1. In the enumerations that follow, 

I = 2, and so the minimum distance of any maximum rank distance (MRD) code, 
i.e. any code that meets the analogue of the Singleton bound, with these parameters 
is 2. As the tables below indicate, there are a number of parameters for which there 
is a self-dual code that is MRD; however, for some parameters, no self-dual codes 
achieve this upper bound. In the latter case, the self-dual codes are all suboptimal, 
and there exist non-self dual codes with the same parameters that are MRdU It is 
not yet clear why/when certain code parameters will not yield any self-dual codes 
that are MRD, but this is an interesting question for future research. 

6. Conclusion and Future Work 

In this paper, we developed a method to enumerate all inequivalent self-dual 
matrix codes of relatively short length over small finite fields.This method is only 
computationally feasible for particularly small parameters because generating a set 
of representatives of double cosets is also computationally intractable for large pa¬ 
rameters. An open question in this area is whether this method may be modified 
to produce an enumeration of self-dual codes over larger fields and/or larger dimen¬ 
sional vector spaces. One possible modification appeals to the existing enumeration 
of block-inequivalent self-dual block codes. It is possible that we may exploit the 
relationship between the block-equivalence and matrix-equivalence maps to gener¬ 
ate representatives of the equivalence classes determined by the intersection of the 
groups of block-equivalence and matrix-equivalence maps using the enumeration of 
block-inequivalent self-dual codes as a starting point. This set of representatives 
would contain all the representatives of the matrix-equivalence classes as well as 
some potentially extraneous codes, but it should be computationally feasible to then 
weed through that list of representatives to find the subset of matrix-inequivalent 
codes. 


tMRD codes can be constructed for any parameters using a modification of a construction from 

[ 11 - 
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Table 1. Enumeration of Linearly Matrix-Inequivalent Self-Dual 
Matrix Codes over F 2 


, Number of Equivalence Class 

q 1 m n ^ 

Equiv. Classes Hepresentatives 


Weight 

Distributions 

2 2 2 4 2 Gi = 

G 2 = 

'1 0 1 o' 
0 10 1 

' 1001 ' 
0 110 

Aq = 1, 2 I 1 = 3 

Ao = 1, Ai = 1, A 2 = 2 

2 2 3 6 5 

‘l 

0 0 10 

o' 


Gi = 

0 

10 0 1 

0 

r- 

II 

T 

T— ( 

II 

0 


0 

0 10 0 

1 



[0 

0 110 

0 ] 


G 2 = 

1 

10 0 0 

0 

II 

T 

m 

II 

T 

t-H 

II 

0 


0 

0 0 0 1 

1 



‘1 

0 0 0 0 

l' 


G 3 = 

0 

10 10 

0 

Ao = 1, Ai = 1, A 2 = 6 


0 

0 10 1 

0 



‘l 

0 0 10 

0 ] 


Gi = 

0 

110 0 

0 

Ao = l,Ai = 5,A2 = 2 


0 

0 0 0 1 

1 



‘l 

0 0 0 0 

1 ] 


G 5 = 

0 

10 0 1 

0 

Ao = 1, Ai = 3, A 2 = 4 


0 

0 110 

0 
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Table 2. Enumeration of Linearly Matrix-Inequivalent Self-Dual 
Matrix Codes over F 3 



1 

m 

n 

Number of 



Equivalence Class 



Weight 

Q 

Equiv. Classes 



Representatives 



Distributions 

3 

2 

2 

4 

1 


Gi 


0 

1 

2 

2 

2 



^0 

= 1, 

A 2=8 

3 

2 

4 

8 

13 

Gi = 

■1 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

2 

0 

2 

0 

2 

0 

2 

0 

1 

O' 

2 

0 

Ao 

= 1, 

A2 = 80 






0 

0 

0 

1 

0 

2 

0 

1 









ri 

0 

1 

1 

2 

0 

2 

21 








G2 = 

0 

0 

1 

0 

0 

1 

1 

1 

1 

0 

2 

0 

2 

0 

1 

2 

Ao 

= 1, 

Ai = 8 , A 2 = 72 






Lo 

0 

0 

0 

1 

2 

2 

oj 









ri 

0 

1 

1 

1 

0 

2 

21 








G3 = 

0 

0 

2 

0 

2 

0 

1 

2 

1 

0 

1 

0 

1 

1 

0 

1 

Ao 

= 1, 

to 

00 

0 






Lo 

0 

2 

1 

1 

1 

2 

2J 









ri 

0 

0 

0 

0 

0 

1 

11 








Gi = 

0 

0 

1 

0 

0 

0 

2 

2 

2 

0 

1 

0 

2 

2 

1 

1 

Ao 

= 1, 

Ti = 4, As = 76 






Lo 

0 

1 

0 

1 

1 

0 

oJ 









ri 

0 

0 

1 

1 

2 

1 

21 








Gs = 

0 

0 

2 

0 

2 

0 

0 

1 

2 

0 

2 

0 

1 

1 

1 

2 

Ao 

= 1, 

to 

00 

0 






Lo 

0 

2 

0 

0 

0 

1 

1. 









ri 

0 

0 

2 

1 

0 

0 

01 








Ge = 

0 

0 

2 

0 

1 

0 

0 

1 

0 

1 

2 

0 

0 

0 

0 

2 

Ao 

= 1, 

As = 80 






Lo 

0 

2 

0 

0 

2 

1 

oJ 









ri 

0 

1 

1 

2 

0 

2 

21 








Gt = 

0 

0 

2 

0 

1 

0 

2 

1 

0 

2 

1 

0 

2 

2 

1 

0 

Ao 

= 1, 

Ai = 16, As = 64 






Lo 

0 

1 

0 

0 

1 

0 

1. 









ri 

0 

1 

1 

1 

0 

2 

21 








Gs = 

0 

0 

1 

0 

2 

1 

1 

1 

0 

1 

1 

0 

2 

0 

1 

0 

Ao 

= 1, 

Ai = 8, As = 72 






Lo 

0 

0 

0 

0 

1 

2 

ij 









ri 

0 

2 

2 

2 

0 

1 

11 








G9 = 

0 

0 

2 

0 

2 

0 

1 

0 

0 

2 

2 

0 

2 

1 

1 

1 

Ao 

= 1, 

Ai = 32, As = 48 






.0 

0 

0 

0 

0 

2 

2 

1 . 









ri 

0 

0 

0 

2 

0 

2 

01 








Cio — 

0 

0 

1 

0 

0 

2 

0 

0 

0 

2 

2 

0 

0 

1 

2 

0 

Ao 

= 1, 

As = 80 






Lo 

0 

0 

1 

0 

2 

0 

1 . 









ri 

0 

2 

1 

2 

0 

2 

21 








Gii — 

0 

0 

1 

0 

0 

2 

1 

1 

1 

0 

2 

0 

2 

0 

1 

2 

Ao 

= 1, 

Ai = 8, As = 72 






Lo 

0 

0 

0 

1 

2 

2 

oJ 









ri 

0 

0 

0 

0 

0 

1 

11 








Gi2 — 

0 

0 

1 

0 

0 

0 

2 

2 

2 

0 

1 

0 

2 

2 

1 

1 

Ao 

= 1, 

Ai = 4, As = 76 






Lo 

0 

2 

0 

1 

1 

0 

oJ 









ri 

0 

1 

1 

1 

0 

1 

11 








Gi3 — 

0 

0 

1 

0 

1 

0 

2 

0 

0 

1 

1 

0 

2 

1 

1 

1 

Ao 

= 1, 

Ai = 20, As = 60 






Lo 

0 

0 

0 

0 

1 

2 

1 . 
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Table 3. Enumeration of Linearly Matrix-Inequivalent Self-Dual 
Matrix Codes over F 4 = F 2 [a] where -|- a -I- 1 = 0 


, Number of 

q 1 m n 

Kquiv. Classes 

Equivalence Class 
Representatives 

Weight 

Distributions 

4 2 2 4 3 = 

■ 1 0 1 01 

_0 1 0 ij 

Ao = l,Ti = 15 

G 2 = 

r 1 0 a 

[ 0 1 a J 

To = l,Ti =3,^2 = 12 

Gs = 

■ 1 0 0 1 1 

_0 1 1 oj 

To = l,Ti =3,T2 = 12 


Table 4. Enumeration of Linearly Matrix-Inequivalent Self-Dual 
Matrix Codes over F 5 



1 

m 

n 

Number of 



Equivalence Class 


Weight 


9 

Equiv. Classes 


Representatives 

Distributions 

5 

2 

2 

4 

2 


G 

- 

;; 

0 

1 

3 

0 

s] 


^0 

= 1. 

Ti = 

24 







G 2 = 

'0 

0 

4 

4 

2 

11 


To 

= 1, 

Ti = 

8, T 2 = 16 

5 

2 

3 

6 

7 

Gi 


'l 

0 

0 

1 

0 

0 

3 

0 

0 

3 

0 ' 

0 

Ao 

= 1, 

Ti = 

124 








0 

0 

1 

0 

0 

3 












'l 

0 

3 

2 

1 

0 ] 










G 2 


0 

1 

4 

1 

1 

1 

Aq 

= 1, 

Ti = 

4, T 2 = 120 








0 

0 

3 

1 

4 

3 












'l 

0 

3 

0 

1 

2 










G 3 


0 

1 

1 

2 

2 

0 

Aq 

= 1. 

Ti = 

4, T 2 = 120 








0 

0 

2 

3 

1 

4 












'l 

0 

2 

4 

2 

0] 










Gi 


0 

4 

4 

1 

4 

1 

To 

= 1. 

Ti = 

28, T 2 = 96 








0 

0 

3 

3 

1 

1 












‘1 

0 

0 

3 

0 

0] 










Gs 


0 

1 

2 

0 

0 

0 

Aq 

= 1, 

Ti = 

12, T 2 = 112 








0 

0 

0 

0 

1 

3 












‘l 

0 

3 

2 

1 

0] 










Ge 


0 

4 

4 

1 

1 

1 

Aq 

= 1, 

Ti = 

28, T 2 = 96 








0 

0 

3 

1 

4 

3 












‘1 

0 

3 

0 

0 

0] 










Gt 


0 

4 

0 

0 

3 

0 

Aq 

= 1, 

Ti = 

44, T 2 = 80 








0 

0 

0 

1 

0 

3 
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Table 5. Enumeration of Linearly Matrix-Inequivalent Self-Dual 
Matrix Codes over F 5 (Continued) 


, Number of 

q 1 m n 

Equiv. C/iasses 



Equivalence Class 



Weight 






Representatives 


Distributions 



5 2 4 8 24 

‘1 

0 

0 

0 

3 

0 

0 

0* 







Gi = 

0 

0 

1 

0 

0 

1 

0 

0 

0 

0 

3 

0 

0 

3 

0 

0 

^0 

= 1, 

Ai = 

624 




Lo 

0 

0 

1 

0 

0 

0 

3J 








ri 

2 

3 

4 

1 

1 

3 

21 







G 2 = 

1 

1 

1 

4 

4 

0 

4 

1 

3 

4 

0 

1 

1 

1 

4 

2 

Aq 

= 1, 

Ai = 

48, 

A 2 

= 576 


.1 

0 

4 

0 

3 

2 

1 

2. 








[2 

4 

1 

3 

3 

2 

1 

41 







G 3 = 

4 

4 

3 

4 

2 

0 

3 

1 

3 

4 

1 

1 

4 

1 

1 

2 

Aq 

= 1, 

.4i = 

48, 

A 2 

= 576 


2 

3 

1 

1 

0 

1 

3 

oj 








[2 

0 

0 

0 

0 

3 

4 

11 







G 4 = 

3 

1 

1 

0 

4 

4 

4 

3 

3 

0 

3 

3 

0 

2 

0 

1 

Aq 

= 1, 

.4i = 

24, 

A 2 

= 600 


Lo 

0 

4 

0 

2 

1 

3 

oJ 








[2 

4 

4 

1 

2 

1 

3 

21 







G 5 = 

1 

1 

3 

3 

1 

2 

2 

0 

1 

0 

3 

1 

0 

1 

0 

2 

Aq 

= 1, 

.4i = 

32, 

A 2 

= 592 


.1 

2 

1 

3 

1 

0 

0 

3J 








ro 

1 

1 

4 

1 

2 

1 

41 







Ge = 

1 

2 

1 

3 

0 

3 

0 

1 

3 

4 

3 

1 

0 

1 

0 

2 

Aq 

= 1, 

.4i = 

64, 

A 2 

= 560 


[3 

3 

2 

2 

2 

1 

3 

oJ 








[3 

4 

1 

3 

1 

1 

3 

21 







G 7 = 

3 

4 

3 

1 

2 

3 

3 

0 

0 

0 

3 

3 

0 

2 

0 

1 

Aq 

= 1, 

.4i = 

24, 

A 2 

= 600 


Lo 

4 

0 

3 

2 

3 

4 

4J 








[4 

4 

0 

2 

0 

1 

3 

21 







Gs = 

4 

2 

1 

0 

0 

1 

0 

0 

0 

2 

1 

2 

4 

4 

1 

4 

Aq 

= 1, 

A 2 = 

624 




L3 

0 

3 

4 

1 

1 

3 

oJ 








ro 

0 

3 

3 

1 

2 

1 

41 







Gg = 

3 

0 

2 

1 

3 

4 

1 

1 

3 

1 

1 

2 

4 

4 

1 

4 

Aq 

= 1 , 

Ai = 

16, 

A 2 

= 608 


L3 

3 

1 

1 

1 

2 

1 

2. 








[4 

1 

3 

1 

2 

4 

2 

31 







G\o — 

3 

0 

4 

3 

4 

0 

3 

3 

1 

1 

3 

1 

0 

1 

0 

2 

Aq 

= 1 , 

Ai = 

24, 

A 2 

= 600 


Ls 

4 

2 

0 

4 

1 

3 

oJ 








ro 

2 

4 

0 

2 

2 

1 

41 







Gii = 

1 

4 

0 

3 

3 

0 

0 

3 

1 

0 

3 

2 

0 

4 

0 

4 

Aq 

= 1 , 

A 2 = 

624 




2 

2 

3 

0 

3 

2 

1 

2. 








[3 

4 

2 

4 

1 

1 

3 

21 







G\2 — 

3 

2 

0 

4 

4 

2 

1 

3 

4 

4 

1 

1 

4 

1 

1 

2 

Aq 

= 1 , 

Ai = 

00 

to 

II 

= 616 


3 

3 

2 

2 

2 

1 

3 

0 
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Table 6 . Enumeration of Linearly Matrix-Inequivalent Self-Dual 
Matrix Codes over F 5 (Continued) 


, Number of 

q 1 m n 

Kquiv. C/iasses 



Equivalence Class 



Weight 




Representatives 


Distributions 

5 2 4 8 continued 

Gi3 — 

■4 

2 

3 

1 

1 

3 

1 

1 

2 

4 

1 

0 

4 

4 

2 

0 

2 

0 

0 

2 

3 

O' 

3 

0 

^0 

= 1, A 2 = 

624 


.2 

4 

3 

1 

4 

2 

1 

2. 





ro 

1 

0 

3 

4 

1 

3 

21 




Gi4 — 

3 

2 

1 

0 

1 

1 

1 

0 

2 

4 

3 

0 

0 

3 

0 

0 

Ao 

= 1, A 2 = 

624 


Lo 

3 

3 

3 

3 

0 

0 

3 J 





ro 

3 

1 

0 

0 

0 

0 

01 




Gi5 — 

4 

3 

0 

3 

0 

1 

2 

4 

4 

3 

4 

2 

3 

4 

2 

4 

Ao 

= 1, Ai = 

144, A 2 = 480 


.2 

2 

4 

1 

1 

0 

0 

3 J 





[2 

4 

1 

3 

2 

2 

1 

41 




Giq — 

3 

4 

4 

4 

2 

2 

1 

3 

1 

2 

3 

1 

0 

1 

0 

2 

Ao 

= 1, A 2 = 

624 


Ls 

4 

0 

3 

0 

4 

2 

1. 





[4 

4 

2 

4 

3 

1 

3 

21 




Gi7 — 

0 

4 

0 

3 

1 

2 

3 

0 

1 

0 

3 

1 

0 

1 

0 

2 

Ao 

= 1, Ai = 

48, A 2 = 576 


Lo 

4 

2 

0 

0 

1 

3 

oj 





ri 

0 

0 

0 

3 

0 

0 

01 




Gie — 

2 

2 

4 

3 

0 

1 

1 

1 

1 

1 

0 

2 

3 

1 

3 

3 

Ao 

= 1, Ai = 

128, A 2 = 496 


0 

0 

0 

1 

0 

0 

0 

3. 





ri 

2 

3 

4 

1 

1 

3 

21 




Gig — 

2 

1 

0 

4 

4 

0 

2 

1 

1 

4 

1 

1 

0 

1 

3 

2 

Ao 

= 1, Ai = 

8, A 2 = 616 


Lo 

1 

4 

2 

0 

1 

2 

3 J 





ro 

1 

1 

2 

2 

0 

3 

11 




G 20 — 

4 

2 

3 

1 

2 

0 

3 

0 

3 

3 

1 

4 

4 

3 

1 

4 

Ao 

= 1, Ai = 

8, A 2 = 616 


L 2 

3 

1 

1 

0 

1 

3 

oJ 





[2 

0 

4 

2 

3 

4 

0 

41 




G 21 — 

0 

4 

0 

0 

1 

2 

3 

2 

1 

1 

3 

2 

0 

0 

0 

1 

Ao 

= 1, Ai = 

16, A 2 = 608 


.2 

1 

0 

0 

4 

2 

2 

4 J 





ri 

4 

2 

0 

4 

3 

0 

31 




G 22 — 

4 

3 

1 

4 

4 

0 

3 

0 

2 

2 

0 

1 

0 

0 

2 

0 

Ao 

= 1, Ai = 

24, A 2 = 600 


.2 

3 

3 

2 

4 

4 

1 

4 J 





ro 

4 

3 

0 

2 

3 

4 

11 




G 23 — 

2 

2 

0 

4 

0 

3 

1 

1 

2 

4 

2 

4 

1 

2 

1 

2 

Ao 

= 1, Ai = 

64, A 2 = 560 


.2 

4 

3 

1 

3 

3 

4 

4 J 





ro 

4 

2 

3 

1 

0 

4 

31 




G 24 — 

4 

4 

3 

0 

4 

4 

0 

1 

1 

3 

1 

3 

4 

0 

1 

2 

Ao 

= 1, Ai = 

144, A 2 = 480 


.1 

0 

0 

1 

3 

0 

0 

3 J 
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Appendix: Details of Implementation of GO„(Fq) in Magma 

Construction of GO„(Fq) for odd q. Magma implements the classical group 
theory definition of the orthogonal group, which preserves a bilinear form that is 
not the dot product; the group that preserves that bilinear form modulo scalars is 
referred to as the conformal orthogonal group GO„(Fq). When n is odd, there is 
precisely one conformal orthogonal group for the specified bilinear form, and this 
group is isomorphic to GO„(Fg) as we define it. When n is even, there are two 
conformal orthogonal groups GO^(F,j) and CO“(F,j), where the plus and minus 
refer to the Witt defect of the underlying bilinear form EH p. 136-141], The 
Witt defect of the dot product determines which of these groups is isomorphic to 
GO„(Fg). Following the background given in [21], we see that 

• if q = 1 (mod 4), then GO„(Fg) is isomorphic to GO(]((Fg), 

• if g = 3 (mod 4) and n = 0 (mod 4), then G02m(Fg) is isomorphic to 
GOL(Fg), and 

• if g = 3 (mod 4) and n = 2 (mod 4), then G02m(Fg) is isomorphic to 
C02-™(Fg). 

Construction of GO„(Fg) for even g. Over characteristic 2, the dot product is 
a singular bilinear form whereas the bilinear form defining COn(Fg) is nonsingular. 
Thus, these groups cannot be isomorphic, and so we resort to a random construction 
of GO„(Fg) in this case. Specifically, we generate a random element A £ GL„(Fg) 
and test whether it satisfies AA^ = A/„ for some A € F*. If A satisfies this condition 
then we include it in the list of generators for GO„(Fg). We continue this process 
until the cardinality of the subgroup generated by this list of generators equals the 
cardinality of GOn(Fg), which can be obtained using formulas for the size of 0 „(F 2 <!) 
from [Tl| and the fact that GO„(F 2 e) is the normalizer of 0„(F2e) in GLji(F 2 e jl. 
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